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In order to explain the anisotropic Rashba-Bychkov effect observed in several metallic surface- 
state systems, we use k-p perturbation theory with a simple group-theoretical analysis and construct 
effective Rashba Hamiltonians for different point groups up to third order in the wavenumber. We 
perform relativistic ah initio calculations for the (V3 x \/3)-R30 o Bi/Ag(lll) surface alloy and from 
the calculated splitting of the band dispersion we find evidence of the predicted third-order terms. 
Furthermore, we derive expressions for the corresponding third-order Rashba parameters to provide 
a simple explanation to the qualitative difference concerning the Rashba-Bychkov splitting of the 
surface states at Au(lll) and Bi/Ag(lll). 

PACS numbers: 71.15.Rf 73.20.At 75.70.Tj 



I. INTRODUCTION 

Since the first experimental verification by LaShell et 
aL— the spin-orbit induced splitting of Shocklcy states on 
metallic surfaces called the Rashba-Bychkov (RB) effect 2 
became into the focus of experimental and theoretical 
research. These investigations range from the prototypi- 
cal L-gap surface states at Au(lll) and Ag(lll)2r— and 
also at Au(110)^ria through Li/W(110) and Li/Mo(110) 
overlayers^i and the Gd(0001) surface^ 2 - to a large num- 
ber of metallic surfaces and surface alloys related to Bi, 
Pb or Sb where the bp and 6p orbitals show a pronounced 
spin-orbit splitting This huge interest is mainly trig- 
gered by potential spintronics applications in relation to 
the Datta-Das transistor,— the spin Hall effect and the 
anomalous Hall effect.— 

While accurate ab initio calculations satisfactorily ac- 
count for most features of the measured dispersion rela- 
tions of metallic surface states, there is an obvious need 
to explain the RB effect in terms of simple models con- 
taining a few, easily identifiable parameters. The sim- 
plest effective Hamiltonian of a two-dimensional electron 
gas, subject to spin-orbit interaction (SOI), includes in 
addition to the kinetic energy, £o + y-^r (k and m* being 
the wavevector and the effective mass of the electrons, 
respectively), a Rashba term,— ^ 

H R (k) = a R (k x a y - k y a x ) , (1) 

where a R is the so-called Rashba parameter and <7, (i = 
x,y,z) denote the Pauli matrices. The corresponding 
eigenvalues, e± (k) = £o + ± a R k {k = |k|) show an 
isotropic splitting for k ^ 0, and, at least for moderate 
values of k, they readily can be fit to most experimental 
and ab initio dispersion relations. 

Although not yet detected experimentally^ a RB split- 
ting that is anisotropic in fc-space is obvious for the sur- 
face states at Au(110). The C2 V point-group symmetry 3 -^ 



of the system not only implies an anisotropy of the effec- 
tive mass, m* ^ m*, but, as discussed in terms of k ■ p 
perturbation theory , 10 ' 35 it leads to a Rashba Hamilto- 
nian containing two independent Rashba parameters, a R 



and cvjj, 



H R (k) = a R k x a v 



(2) 



Fully relativistic ab initio calculations confirmed the exis- 
tence of the anisotropic RB splitting at Au(110),— match- 
ing with a high accuracy to the eigenvalues of the effective 
Hamiltonian in Eq. ([2]) J£ 

Even in case of high symmetry surfaces, i.e., having a 
point-group of C 3v or d v , several studies^ 3 -^^^^^ 
called the attention to an anisotropic RB splitting. In 
Ref. [13 the anisotropic RB effect at Bi/Ag(lll) and 
Pb/Ag(lll) surfaces was reproduced by using a nearly- 
free electron model and explained due to in-plane struc- 
tural inversion asymmetry. From the group-theoretical 
analysis in Ref. |3fj| it is, however, clear that under C% v 
and C± v point-group symmetry an effective 2x2 Hamil- 
tonian that is linear in the components of k must be of 
the form of Eq. ([I]), hence, it can not explain the observed 
anisotropy of the RB splitting. Thus, we conclude that in 
these systems the anisotropic RB effect can be described 
by a Hamiltonian containing at least third-order polyno- 
mials of k x and k y . It should be noted that the second- 
order terms are related to the kinetic energy (effective 
mass terms) that are irrelevant to the RB splitting. 

To construct Rashba Hamiltonians up to third order in 
k, in the present work we use k-p perturbation theory and 
group-theoretical methods different from Ref. [35l Our 
analysis of the effective Hamiltonian is closely related to 
that of Ref. [HI, where, for the case of Cs v symmetry, the 
correct form of HQs) is derived up to third order in k 
and the corresponding band dispersion was used to ex- 
plain the hexagonal warping of the surface states' Fermi 
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contour observed experimentally in the topological insu- 
lator Bi 2 Te 3 .^ 

We also perform relativistic ab initio calculation for 
the Bi/Ag(lll) ordered alloy in (-y/3 x y/S)B30° super- 
structure and confirm that for higher values of k, but still 
in the measured range, third-order terms of the Rashba 
Hamiltonian are needed to reproduce the RB splitting 
and that these terms are of the predicted form. More- 
over, using explicit expressions of the third-order Rashba 
parameters within k ■ p perturbation theory and calcu- 
lated spectral densities at the Brillouin zone center we 
are able to give a simple explanation why the Au(lll) 
and the Bi/Ag(lll) surface states exhibit isotropic and 
anisotropic RB splitting, respectively. 



II. PERTURBATION THEORY AND 
SYMMETRY ANALYSIS 



Let Q be a high symmetry point of the surface Bril- 
louin zone (SBZ), for which a pair of spin-degenerate 
cigenstatcs exists on a nonmagnetic surface. Due to time 
reversal symmetry, this is always the case if Q = — Q + K 
is satisfied where K is a two-dimensional (2D) reciprocal- 
lattice vector. Such points are the center of the SBZ (r) 
and some special points at the boundary of SBZ such as 
the X, Y and S points for a primitive rectangular lat- 
tice, the M and X points for a square lattice and the M 
point for a hexagonal lattice. As what follows, our in- 
vestigations will concern solely this case termed as the 
proper Rashba effect £ As pointed out in Ref. [tl due to 
double-group symmetry there can happen degeneracy at 
points of the SBZ that doesn't meet the above condition, 
like the K point for a hexagonal lattice (improper Rashba 
effect). 

To describe the surface band around Q it is worth to 
label the corresponding Bloch states by the wavenumbcr 
with respect to Q, V'Q+k, and introduce a new wavefunc- 
tion, 0k as 



and 



= iQT 



(3) 
0k (r), 



V Q+k (r) = e* kr k (r) , 

with the boundary condition, 0k (r + T) 
where T is a 2D real-lattice vector. Considering the 

2 

Hamilton operator, % = j— + V + Hso, with the crystal 
potential, V, and Hso denoting the spin-orbit interac- 
tion, 



Hso 



h 



(W x p) • er 



4m 2 c 2 

the wavefunctions 0k satisfy the eigenvalue equation, 



(4) 



(Ho(k)+H so (k))0 k = e(k)0k, 



with 



(fik + pf 
2m 



+ V, 



(5) 



(6) 



^so(k) = 



4m 2 c 



(W x (hk + p)) ■ er 



(7) 



Following the recipe used in Ref. [T(| in the first step 
we look for the solution of the Schrodinger equation, 



Ho(k)0 k =e o (k)0k 



k ! 



(8) 



which can be elucidated, e.g. in terms of k ■ p perturba- 
tion theory. Although such a calculation provides with 
a deeper insight into the problem^ in this section we 
just make use of the symmetry properties of the solu- 
tions, £o(k) and k . First we note that since 'Ho(k) 
is independent of the spin, the solutions of Eq. ([SJ re- 
main degenerate in spin-space. Time reversal symmetry, 
TH (k)T- 1 = Ho(-k) with Tip = tp*, then immediately 
implies 



£ (-k) = e (k) , 



(9) 
(10) 



where the phase of k has been fixed without loss of gen- 
erality. Clearly from Eq. ((9]), a polynomial form of £o(k) 
contains just even powers: the first non-trivial (second- 
order) terms are obviously related to the effective masses. 

We can draw further relations from point-group sym- 
metry. Let Gq be the small group of Q, i.e. gQ = 
Q + K for any g £ Gq and K denoting an appropri- 
ate reciprocal-lattice vector. Using the standard defini- 
tion for the action of a symmetry operation, (g o /)(r) = 
/(5 _1 r), from the symmetry of the Hamilton operator, 
g o %o(k) = Ho(gh), one easily can derive 



e (gk) = £ (k) , 

g°4>l 



(11) 

(12) 



In the second step, using "Hso(k) as perturbation and 
^Xs with Xs being spinor eigenfunctions (s = ±i) as 
unperturbed wavefunctions, first-order degenerate per- 
turbation theory is applied. The Rashba Hamiltonian, 
Hft(k), is defined as the corresponding 2x2 matrix, 



where 



<(k) = (0 k | 



H R (k) = a(k) • a , (13) 



(Wx (ftk + p))|0°) . (14) 



Our present goal is to derive the polynomial form of 
ct(k). To this end we note two symmetry properties that 
can be obtained from Eqs. (fT0|) and (|T2|) : 



(15) 



a(-k) = -a(k) 



stating that cti(k x , k y ) can be expanded in terms of poly- 
nomials of odd power and 



a(g k) = det(g) ga(k) 



(16) 
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where det(g) = 1 for proper rotations and det(g) = — 1 
for improper rotations. Eq. (fTB) is then used to set 
up linear equations for the coefficients, c\, of the n th - 
order polynomials of ai(k x ,k y ) = Ylt=i n c Wx^y~ l 
(i = x, y, z). Solving this set of linear equations serves to 
search for the vanishing coefficients, in principle, for any 
power n, hence, to determine the form of Hn(k). 

Another systematic way to obtain Hn(k) relies on the 
observation that Eq. (fT!)|) can be used to formulate the 
invariance of the Rashba Hamiltonian as 

H R (k) = ot(gk)-(det(g)ga) , (17) 

also implying that <x transforms as an axial vector. Sort- 



ing out the components of k and er according to irre- 
ducible representations of Gq, their direct products can 
again be decomposed into irreducible representations. 
Eq. (|f 7p states that only the total symmetric irreducible 
representations from this decomposition can contribute 
to -fffj(k). From the corresponding tables of the point 
groups^ one can easily construct the possible terms en- 
tering iJfj(k) according to increasing powers of k x and 
k y . In the Appendix this procedure is illustrated for the 
simple case of point group Ci v . As one of the main result 
of this work, in Table U we list the possible terms up to 
third order in k that can enter H R (k) for different point 
groups relevant to surfaces of crystals. 
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k x &X H~~ kyCJy , 

^x^y kyO~x ; 
k x ky(j x k x ky(jy , 

kxkyCJ x ~\~ k x kyO y 


k x Oy i 

k x ky(J x ) 

kxkyCy 3 

ky(J X 


{k x ~\~ kxky) (jy 
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TABLE I. Possible terms of the Rashba Hamiltonian for different point groups (first row) containing first-order (second row) 
and third-order (third row) polynomials of k x and k y . 



Finally in this section, we comment on the method 
used in Ref.Uft In this work a Hamiltonian including SOI 
but excluding all fc-dependent terms was considered as 
the unperturbed system and the twofold degenerate solu- 
tions, (f>\ and 02, corresponding to the wavenumbcr Q as 
the unperturbed solutions. The perturbation was there- 
fore taken as H'(k) = ^kp + -^r^ (VV x k) • a and, 
similar to our strategy, first-order degenerate perturba- 
tion theory was applied. The form of the effective Rashba 
Hamiltonian, H^(k) = (fa \ H'(k) \ fa) (i,j = 1,2), is 
then determined via the invariance conditions, 

H'(gk) = D(g) H'(k)D(g)-\ (18) 

where D (g) is a 2 x 2 unitary double-point group rep- 
resentation of g. In case of Abclian point groups (Chx, 
C2, C3 and C4), the degenerate states form time- reversed 
pairs and D(g) can simply be set up from the characters 
of the corresponding one-dimensional irreducible repre- 
sentations. Following from the definition of H'(k), in 
Ref. HH the first-order Rashba Hamiltonians were ob- 
tained for the groups Chx, C21,, C^ v and C± v . It is, how- 



ever, straightforward to show that, when applied to the 
Hamiltonian (fT3]) , Eq. (fT8|) is equivalent with condition 
(fl~7|) .— hence, using double-group representations leads 
to the same results as listed in Tabic Q] 



III. THIRD-ORDER RASHBA SPLITTING AT 
BI/AG(111) 

By using the relativistic Screened Korringa-Kohn- 
Rostokcr (KKR) method^ we performed calculations for 
the (V3 x ^)i?30° ordered surface alloy Bi/Ag(lll) to 
obtain a quantitative verification of our prediction of 
a third-order Rashba Hamiltonian. A 2D lattice con- 
stant of 2.892 A related to fee Ag bulk and, accord- 
ing to geometry optimization we performed in terms of 
the VASP method^ and also in agreement to previ- 
ous LAPW calculations^ an outward buckling of 36 % 
(0.85 A) for the Bi atoms were considered. The local 
spin-density approximation as parametrized by Vosko et 
al^ was applied, the effective potentials and fields were 
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treated within the atomic sphere approximation with an 
angular momentum cut-off of £ m ax = 2. The energy in- 
tegrations were performed by sampling 12 points on a 
semi-circular path in the upper complex semi-plane and 
for the necessary ^-integrations we selected 36 Appoints 
in the irreducible segment of the surface Brillouin zone. 

The calculated dispersion relation of the Bi surface 
states below the Fermi level is shown in Fig. [T] along 
the r — M direction of the SBZ. The maxima of the 
Rashba-split Bi sp z band are shifted from the T point by 
fco = 0.1 A -1 and, using a parabolic fit around the max- 
ima, we obtained an effective mass of m* = — 0.36m e . 
These values are in good agreement with experimental 
data, fco = 0.13 A -1 and m* = — 0.35 m e ^ It should, 
however, be mentioned that the calculated surface bands 
are shifted downwards by about 0.5 eV as compared to 
experiment, most probably, due to the angular momen- 
tum cutoff and to the atomic sphere approximation used 
in the calculations. The Bi p x p y surface bands, shifted 
upwards due to crystal field splitting and to spin-orbit 
coupling, can also be clearly seen in Fig. [T] 
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FIG. 1. Calculated dispersion relations of the occupied sur- 
face states of Bi/Ag(lll) along the V — M direction of the 
SBZ. 



_ InFig. |we plotted Ae(k) 2 along the T - K and the 
T — M directions, together with different fitting functions 
related to Eq. (f2U|) . It can be seen that a parabolic fit 
(dots), otik 2 with a.\ = 1.74 eV A, applies well to the two 
curves only for about k < 0.07 A -1 . Up to k ~ 0.13 A -1 
the two curves still coincide, however, the isotropic third- 
order contribution is needed for a good fit: here we used 



fc 3 ) 



with the same value for 

3 



a fitting function (aifc + a 3 

For wavenumbers 
k > 0.13 A -1 , the anisotropy of the RB splitting be- 
comes apparent: along T — M ((f) = ir/2) the previous fit 
applies, while along T — K (0 = 0) the fitting function 

had to be extended by the anisotropic term, (a§) fc 6 

with a§ = 9.4 eV A 3 . Our numerical results thus clearly 
support the appearance of third-order terms for the sur- 
face states of Bi/Ag(lll) consistent with the functional 
form as derived from group-theoretical methods. 
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FIG. 2. Square of the calculated splitting, Ae(k) = (e+(k) — 
s_(k))/2, of the occupied surface states of Bi/Ag(lll). 
Squares: T — K direction, circles: T — M direction, see the 
sketch of the SBZ in the inset. Dotted and solid lines display 
first-order and third-order fits as described in the text. 



In case of C^ v symmetry the effective Rashba Hamil- 
tonian can be written up to third order in k as (see Table 

ID, 

HrQc) = (a t + a\ fc 2 ) (h x a y - k y a x ) + a 2 fc 3 cos 3cp a z , 

(19) 

with the polar coordinate <p — &iccos(k x / k) . Note that 
the x axis was chosen along the T — K direction of the 
SBZ. Obviously, there are two kinds of third-order con- 
tributions to the Hamiltonian (|19[) : an isotropic one with 
coefficient 0:3 and an anisotropic one with the coeffi- 
cient a§. The square of the splitting of the eigenvalues, 
Ae(k) = (e+(k) — e_(k))/2, can then be expressed as 

Ae(k) 2 = (aifc + a\ fc 3 ) 2 + (a 2 ) 2 fc 6 cos 2 3cp . (20) 



IV. COMPARISON OF THE RASHBA EFFECT 
AT AU(lll) AND BI/AG(111) 

It is well-known from experiments and ab inito 
calculations that the Au(lll) i-gap surface states 
show a highly isotropic (first-order) Rashba splitting. 
Since both systems, Au(lll) and Bi/Ag(lll), exhibit 
Cs v symmetry, the question naturally arises why there is 
a remarkable difference concerning third-order RB split- 
ting. In order to find, at least, a qualitative understand- 
ing of the problem we extended the fc • p perturbation 
calculations presented in Ref. [13 for the case of Gi v sym- 
metry to C% v symmetry and found the following expres- 
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sions for the third-order coefficients in Eq. (JTHJ) , 



4m 4 c 



and 



E 



E 



(21) 



{UP X \^i)^t\d X V\^ n ){^ n \pM 

(e ~ e^)(e - e£) 



(22) 

From the above formulas it turns out that third-order 
corrections to the effective Hamiltonian arise from an 
admixture between the surface state, 4>o of sp z orbital 
character at energy eq, and those corresponding to the 
two-dimensional irreducible representation, E, (f>^ with 
p x ± ip y character, at energy e^. Note that all these 

states are eigenstates of the Hamiltonian, % = ^— + V, 
at the center of the surface band Q. It is remarkable that, 
similar to the isotropic first-order Rashba parameter, the 
strength of the isotropic contribution, a\, depends on 
the partial derivative of the crystal potential normal to 
the surface, d z V, while the coefficient for the anisotropic 
term, a| , is related to the in-plane gradient of the poten- 
tial, d x V. 

In Fig. |3] we plotted the scalar-relativistic, orbital pro- 
jected partial densities of states (Bloch spectral func- 
tions) at the r point in an appropriate energy window 
around the surface states of Au(lll) and Bi/Ag(lll). 
From the upper panel it can be seen that in case of 
Au(lll) the edge of the bulk subband of E symmetry 
closest to the surface state is by Ae = 1.77 eV below eq. 
The situation is entirely different for Bi/Ag(lll), since 
the Bi p x py states of E symmetry are separated from the 
sp z state by only Ae = 0.27 eV. Since these are the char- 
acteristic energy differences that enter the denominators 
in Eqs. (|2"Tj) and (f2"2"j) . a difference of at least two orders in 
magnitude can indeed be estimated concerning the third- 
order RB effect. Most probably, the actual values of the 
matrixelemcnts of p x and d XtZ V even further strengthen 
this difference. 



V. CONCLUSIONS 

Based on k ■ p perturbation theory including spin-orbit 
interaction we gave a suitable definition to an effective 
Hamiltonian, Eqs. (fT3")) and (O, describing the Rashba- 
Bychkov splitting on metallic surfaces. Due to time rever- 
sal and point-group symmetry we showed how to obtain 
the most general forms for the effective Hamiltonian, and 
derived them up to third order in k for point groups com- 
patible with surfaces of real crystals. Since the effective 
Hamiltonan (|I3[) applies to a couple of non-interacting 
two-band models, the expressions listed in Table U can 
be used in quite a general sense. 

Using the rclativistic Screened Korringa-Kohn- 
Rostokcr method, we demonstrated that the Rashba 
splitting of the Bi sp z surface band of the ordered surface 
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FIG. 3. Calculated orbital projected partial densities of states 
at the r point around the surface states of Au(lll) (upper 
panel) and Bi/Ag(lll) (lower panel). 



alloy Bi/Ag(fff) can not be satisfactorily described 
in terms of a first-order isotropic Rashba Hamiltonian. 
Moreover, we showed that the strong third-order contri- 
bution is subject to an anisotropy consistent with the 
dispersion relation deduced from our symmetry analysis. 

We also derived explicit formulas for the third-order 
anisotropy parameters and established that the isotropic 
and anisotropic contributions are related to the normal- 
to-plane and the in-plane gradients of the crystal po- 
tential, respectively. Comparing the energy separation 
of relevant orbital projected bands for Au(fff) and 
Bi/Ag(fff), the derived expressions were useful to give 
a qualitative understanding of the different nature of 
Rashba-Bychkov splitting in these two systems. 
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Appendix: Rashba Hamiltonians for point group C-2v 

By using the direct products of irreducible represen- 
tations, in this Appendix we give an example for the 
polynomial forms of a 2 x 2 effective Hamiltonian for the 
point-group C 2v ■ Let us denote the elements of the group 
by E : {x,y,z}, C 2 : {-x,-y,-z}, S x : {~x,y,z} and 
S y : {x, — y, z}. The group has four one-dimensional ir- 
reducible representations: A±, A 2 , B±, B 2 with the char- 
acter tabled 





E 
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Sy 
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-1 
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-1 


-1 
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Making use that the k and <x transform as polar and 
axial vectors, respectively, a comparison with the char- 
acter table lets us to sort out the components of this vec- 
tors according to irreducible representations: B\ : k x ,a y , 
B 2 : k y , a x and A 2 : a z . 



From the tabic of direct products, 




Ax A 2 Bi B 2 


Ai 


A x A 2 B x B 2 


A 2 


A 1 B 2 Bi 


B! 


A-l A 2 


B 2 


Ay 



it is easy to find that the only combinations that are first 
order in k x and k y and correspond to the A\ irreducible 
representations are k x a y and k v <r x , therefore, the first- 
order Rashba Hamiltonian can be written in the form of 
Eq. ©. 

The second-order polynomials of k x and k y can be 
sorted according to irreducible representations as follows: 
Ay : k 2 ,, ky and A 2 : k x k y . It should be noted that this 

implies the form of Tp-^-kZ + Tp^k 2 for the effective mass 

^ 2m* x 2m* V 

term. The third-order polynomials of k x and k y can then 
be classified as Bi : kl,k x k 2 and B 2 : ky,k y k 2 . Tak- 
ing direct products with cr^ of A\ symmetry leads to the 
possible third-order contributions to the Rashba Hamil- 
tonian: k^cTy, k 2 k y cr x , k x kyO y and ky(T x , i.e. 

fl|.(k) = a\ k'la y + a% k 2 x k y a x + a\ k x k 2 a y + a\ k z u x . 

(A.l) 
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